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The question of whether it is possible to suppress ray oscillations in light 
waveguides is important for the design of light communications systems. 
With the help of Liouville's theorem of statistical mechanics it is shown that 
it is impossible to reduce simultaneously the amplitudes and the angles of 
ray oscillations if the ray originates in and returns to a region of low index 
of refracti(M. A reduction of both ray amplitudes and angles can be achieved 
only if the ray moves from a region of low to one of high index of refraction. 

Liouville's theorem is used to derive a condition relating the output posi- 
tion and slope of a ray which traverses an optical transformer to its input 
position and slope. With p, , x, denoting the canonically conjugate variables 
of the output ray and Pi , xi those of the input ray, the condition derived 
from Liouville's theorem states that the Jacobian of the transformation is one. 

d(Vi , X.-) _ ^ 
d{pi,Xi) 

1. INTRODUCTION 

Light transmission systems can be built in various ways. A continuous 
dielectric medium of rotational symmetry with an index of refraction 
which depends on the distance r from the optical axis 

n = n{r) 

is capable of guiding light rays if n(r) decreases monotonically with in- 
creasing r. Another example is the beam- waveguide consisting of a 
series of lenses which refocuses the Hght beam periodically counteracting 
diffraction. 

Both of these examples have one point in common — a ray which i.s 
launched off-axis into the waveguide follows an oscillatory trajectory. 
However, even if a light ray travels on-axis it will he forced into an 
oscillatory trajectory by any imperfection of the guidance medium.' To 
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keep the ray amplitudes small requires a very high precision of ahga- 
ment which might be hard to obtain for long waveguides. 

It seemed natural, therefore, to consider means of suppressing these 
ray oscillations, and if all such efforts fail, to ask for a general physical 
principle which says that such ray oscillation suppressors are impossible. 

The search for such a general principle is even more important as it is 
easy to construct models of beam waveguides which violate physical 
principles in subtle ways thus seeming to lead to ray oscillation suppres- 
sors. One such system is shown in Fig. 1. Assume that we deform thin 
lenses as indicated in the figure and assume further that these lenses 
behave just like plane thin lenses in that they break each ray by an 
amount ^„ which depends only on the radius r„ of the ray but not on 
the input angle. 

tan/3„ = -r„//. 

Making the paraxial approximation, which means replacing tan (3„ by 
0„ and tan a^ by «„, we obtain the ray equation 

r„+i = r„ -f a„(2„+i — Zn) ~ '' (la) 

«„+i = «„ ^. (lb) 

If the lenses are warped to form parabolas, 

z„+, - z^= d + 6(r„2 - Tn+i'). (2) 

Equations (la) and (lb) together with (2) describe rays which, if they 
travel from the left to the right in Fig. 1, exhibit decreasing amplitudes. 
In fact, if one allows each ray to travel a sufficient distance they approach 
the axis arbitrarily closely. 

It appears that we have invented a ray oscillation suppressor. 

The object of this paper is to prove that such a device is impossible. 
So the question arises : What went wrong with the argument presented 
above? A closer examination shows that the assumption that j3„ is 




LENS: n n+r 

Fig. 1 — Beam-waveguide composed of warped, thin lenses. 
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independent of a„ violates Liouville's theorem. We will return to this 
question later. 

The general proof of the impossibility of constructing a ray oscilla- 
tion suppressor was suggested by J. R. Pierce. 

II. PROOF OF THE IMPOSSIBIUTY OF A RAY OSCILLATION SUPPRESSOR 

The proof is based on Liouville's theorem.^ It refers to the representa- 
tion of physical systems in phase space. Phase space is the space of the 
canonically conjugate variables 5, and p, describing the system. Each 
system is represented by one point in phase space. Many identical sys- 
tems which happen to be in different states described by different values 
of their coordinator g^ and Pi can be described by the density of their 
representation points in phase space. Liouville's theorem states that the 
density of any given configuration of points in phase space is constant 
if the systems under consideration obey the canonical differential equa- 
tions 

^" = ^ ^' = — ^ (3) 

dt dpi dt dqi ' 

H is the Hamiltonian function describing the system. Another version 
of Liouville's theorem states that the volume containing a constant 
number of points in phase space remains constant in time. 

For Liouville's theorem to be applicable to light rays we have only to 
show that light rays can be described by equations of the form (3). The 
derivation of the Hamiltonian equations of geometric optics can be 
found in Ref. 5. The derivations are sketched here for the sake of con- 
venience. 

To show this we start with Fermat's principle which states that a 
light ray connecting two arbitrary points Pi and Pz in a medium of index 
of refraction 

n = n(x,y,z) (5) 

follows a path such that 

J = ~ I nd s = extremum. (6) 

c Jpi 

Here, c is the velocity of light in vacuum and s is the path length meas- 
ured along the ray trajectory. Introducing coordinates x,y,z, we can 
rewrite (G) as 

J = - n Vl + x"^ + y"^ ^2 ^ extremum (7) 

c Jpi 
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with 

t dx , , dy ,„v 

X = ^- and y = -j-. \o) 

dz dz 

Equation (7) is analogous to Hamilton's principle of least action with 
the Lagrangian 

L = n Vl + x'' + y'^ (9) 

and the time t being replaced by the z-coordinate. 

Once the Lagrangian of a system is known the moments p^ and py 
canonically conjugate to the x and y coordinates are defined by 

(10a) 

(10b) 
VI -t- 3;' -I- y 

and the Hamiltonian function by 

H = pxX + pyy — L = — "Vn — p/ — py. (U) 

The variational problem (7) is solved by the equations^ 

p. = - ^ p. = - a^- (12b) 

Equations (12a) and (12b) are analogous to (3) which shows that the 
ray description can be given in terms of canonical differential equations. 
The equations of (12a) are satisfied identically while the equations of 
(12b) lead to the well-known ray equations 





dL 
dx' 

dL 

dy' 


X 


P' 


'"Vi-tx'^ + y" 

y' 


py 


'Vl+x'^ + J/'^ 



1 d_ ( x' \ ^ dn 

Vl + x'2 + y'^ dz V Vl + x'^- + y'y dx 

1 d / y' \ ^ dn 

Vl + x'^ + y" dz V Vl + x'^ + y'y dy- 

Introducing 



(13a) 
(13b) 



^' - Vl + x"^ + y'\ (14) 

dz 

(13a) and (13b) can bo written in the more famiUar form* 
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(Is \ dsj 



ds 



dn 
dn 



(15a) 
{15b) 



The preceding discussion of ray dynamics was sketched only to prove 
that Liouville's theorem applies to light rays. 

Now, we are finally in a position to prove the impossibility of a ray 
oscillation suppressor. To simplify the discussion let us limit the problem 
to two dimensions, x and z. Assume that z is the axis of the system. The 
phase space is now two dimensional and is spanned by the coordinates 
X and px . Let us further assume that we consider an ensemble of rays 
whose initial conditions are such that the representation points of all 
these rays fill a square area centered around the origion of phase space 
as shown in Fig. 2. Each ray represeiited in this area has a certain dis- 
tance X from the optical axis z and a certain slope given by (10) 



X ^ 



■y/n^ — pj^' 



(16) 



If an oscillation suppressor were possible we would require that all the 
rays initially contained in the square of phase space of Fig. 2 would 
approach the z-axis more closely. In addition, we would require that 
the angles between the rays and the s-axis don't increase or perhaps 
even decrease. If we look at the rays initially and finally in a region of 
constant index of refraction n for example in vacuum, n = 1, we would 
find that the square of Fig. 2 has deformed either into the rectangle, if 
the angles don't shi'ink, or into the smaller square, if the angles as well 
as the amplitudes shrink, as indicated in Fig. 2 by dotted lines. In either 



Fig. 2 — Volume in phase space occupied by light ray representation points. 
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case, we find that the area (volume in two dimensions) of phase space 
occupied by the points representing the initial ray positions has de- 
creased. However, Liouville's theorem says that this is impossible so 
that we see that a ray oscillation suppressor is impossible. We can trade 
off amplitude at the expense of spread in -p,, direction. In this case, either 
the tangent of the ray angles x' or the index of refraction has to increase. 
It is even possible to decrease both the ray amplitudes and angles by 
increasing n along the 2-axis. However, the area in phase space has to 
stay constant. The initial square has to deform into a rectangle of equal 
area which stretches along the px axis. After some distance we have 
reached a region of high index of refraction and find both the amplitude 
and the ray angles decreased (but not the -px values which have in- 
creased). For most appUcations to be able to make use of the effect, we 
have to leave the high index medium. But as soon as n drops to a low 
value the angles have to increase to keep the spread in -pi constant and 
again we have traded a decrease in the ray amplitudes for an increase 
in the ray angles. The ray position in most optical systems will eventually 
spread far apart if we allow the rays to travel far enough. This means 
that the volume in phase space, though its volume content remains 
constant, assumes a "filamentous" appearance and extends to many 
different parts of phase space.^ 

in. A BASIC RELATION FOR OPTICAL THANSFORMEHS 

Liouville's theorem allows one to formulate a theorem which all rays 
passing through an optical device (optical transformer) have to obey. 

Let us assume we have an arbitrary optical transformer with input 
rays whose positions and slopes are described by the canonieally con- 
jugate variables g,- , p< and corresponding output ray with variables 
q,- , p, . The output variables are related to the input variables by 

q> = qi(2i , p.) 
P.- = Pi(s.- , Vd- 

The input rays may occupy a volume dV ^ dqidq^idpidpz in phase space. 
This volume deforms, as the rays propagate, to dV. Liouville's theorem 
states that these volumes are identical; 

dV = dV. (17) 

The volume dV on the right hand side of (17) can be rewritten as 

dV = l^^i^ , dq, dq, dp, dp, (18a) 
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or 



dV = ^P^l dV. 



(18b) 



Wo conclude from {17) and (18) that the Jacobian must be equal to 
unity 

d(qi, p,) 



d{qi, V<. 



= 1. 



(19) 



Equation (19) is stated in Ref. without proof. 

The derivation of (19) is based on the fact that the ray trajectory can 
be described by the differential equations of (13). However, there may 
be discontinuities in the index of refraction, n, where the ray equations 
can not be applied. But it is well known that rays traverse discontin- 
uities of the index of refraction. The ray trajectory is unaltered if the 
discontinuity is replaced by a rapidly changing but continous transi- 
tion of n. In this way we assure that the ray equations hold everywhere 
and that (19) is applicable even in that case. 

Limiting the problem to two dimensions we can write (19) as 



ap ax 

dp d-v 



_ dpdx _ 
dx dp 



(20) 



Equation (20) allows us to derive an interesting relation between the 
input and output angles of rays passing through an intinitesimally thin 
optical transformer (Fig. 3). If the thickness of the optical transformer 
shr inks to zero we have x = x and consequently, 




OPTICAL -— 
TRANSFORMER 




Fig. 3 — Illustrations of a thin optical transformer. 
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dp dx 



SO that (20) reduces to 

dp 



^ = 1 



whose solution is 

V = p-^f(x). (21) 

With the help of (10) we see that if 

dx , dx. , 

= tan a -j- = tan a, 



dz dz 



it follows that 



p = n sin a p = n sm a 

so that (21) can be written as 

71 sin a = n sin a + (n sin a)a^a ■ (22) 

This is a fundamental relation which all rays passing through thin lenses 
or any other thin optical device have to obey. 

If both a « 1 and a « 1 and n = n = 1, (22) simplifies 

^ a - a =^ (5)_o . (23) 

This is the relation which is used to describe the change of ray angles 
passing through a thin lens. Equation (22) shows that this thin lens 
relation holds approximately for rays which impinge nearly perpendicu- 
lar to the lens. If the rays make large angles with respect to the direction 
normal to the lens surface (23) has to be replaced by (22) . This explains 
the error which was made in deriving (1). If this equation is corrected 
by using (22) rather than (16) the ray oscillation suppressing quality 
of the warped thin lenses disappears. 

The general expressions (19), (20), or (22) can be used to check the 
physical realizability of optical models. 
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